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ABSTRACT: We demonstrate a close connection between the decoherent histories 
(DH) approach to quantum mechanics and the quantum state diffusion (QSD) picture, for 
open quantum systems described by a master equation of Lindblad form. The (physically 
unique) set of variables that localize in the QSD picture also define an approximately 
decoherent set of histories in the DH approach. The degree of localization is related to the 
degree of decoherence, and the probabilities for histories prescribed by each approach are 
essentially the same. 



The last decade has witnessed considerable interest in the foundations of quantum 
mechanics. Many reasons for this may be found: the long-felt dissatisfaction with the 
Copenhagen interpretation; certain experimental developments and the Copenhagen in- 
terpretation's inability to supply a useful qualitative account of them; the special needs of 
quantum mechanics applied to the entire universe (quantum cosmology); and the general 
desire to possess a deeper understanding of quantum theory. Partially for these reasons, 
modifications and generalizations of both the mathematics and interpretation of quantum 
theory have been sought. 

This letter is concerned with demonstrating the connections between two recently 
developed alternative approaches to quantum theory, each of which was proposed, inde- 
pendently, with the aim of shedding light on some of the difficulties outlined above. The 
approaches we shall compare are the decoherent (or "consistent") histories approach [1,2] 
and the quantum state diffusion picture [3,4,5] (see Ref.[6] for an early guess at the rela- 
tion between these two approaches). In both of these approaches, the basic mathematical 
formalism of quantum mechanics is left untouched, but new insight into its interpretation 
is obtained by focusing on different types of mathematical objects. 

We will be concerned with a quantum system consisting of a subsystem coupled to its 
environment. The subsystem is then frequently referred to as an open quantum system, 
and we shall do so here. Conventionally, an open system is described by a reduced density 
operator p, evolving according to a master equation, derived by tracing over the environ- 
ment. Under the assumption that the evolution is Markovian, the master equation takes 
the Lindblad form [7], 



Here, H is the Hamiltonian of the open system in the absence of the environment (some- 
times modified by terms depending on the Lj) and the n operators Lj model the ef- 
fects of the environment. This equation is frequently used in quantum optics [8] and 
in studies of decoherence [9,10]. For example, in the much-studied quantum Brown- 
ian motion model [11,12,13], the master equation is (1) with a single Lindblad operator 
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L = (2D)-2(x + 2^Dp), with D = h 2 /(8m~fkT) (where 7 is the dissipation and T is the 
temperature of the environment), and H = Hg + ^{x,p}, where Hg is the distinguished 
subsystem Hamiltonian, in the absence of the environment. 

Density operators satisfying (1) give statistical predictions in full agreement with exper- 
iment in a wide variety of situations. However, they do not give a picture of the behaviour 
of an individual system, but only of ensembles. The quantum state diffusion approach 
avoids this shortcoming [3,4,5,14]. It originated from considerations of the quantum mea- 
surement problem [15,16,17,18], and the desiderata of describing individual experimental 
outcomes and putting physical intuition into the equations, as urged by Bell [19]. It was 
also motivated by its computational advantage and insight in treating practical problems 
in open systems [20]. It consists of an "unravelling" of the evolution of p under (1). 
This involves regarding p as a mean over a distribution of pure state density operators, 
p = M\ip){ip\, where M denotes the mean (defined below), with the pure states evolving 
according to the non-linear stochastic Langevin-Ito equation, 
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+ J2(L j -(L j ))\i>)dt j (t) (2) 
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for the normalized state vector Here, the d^j are independent complex differential 
random variables representing a complex Wiener process. Their linear and quadratic 
means are, M[dfjdf£] = S jk dt, M[d^jd^ k ] = and M[d£j] = 0. 

The quantum state diffusion picture described by the Ito equation (2) is mathematically 
equivalent to the Lindblad equation (1). However, its appeal lies in the fact that the 
solutions to the Ito equation appear to correspond rather well to individual experimental 
runs, and thus provide considerable insight into the behaviour of individual processes 
and systems. Solutions to the Ito equation commonly have the property of localization 
- the dispersion of certain operators tends to decrease as time evolves. This has been 
demonstrated by numerical solutions [5], analytic solutions in special cases [13,18,21], and 
some general theorems [4,14]. The method has also been successfully used to analyze 
quantum jump experiments [22]. 
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Our first results, required for the comparison with the decoherent histories approach 
below, are explicit representations of the solutions to the Lindblad equation (1) and the 
Ito equation (2). To solve the Lindblad equation, consider the case of a single Lindblad 
operator L = Lr + iLj, where Lr, Lj are hermitian. Divide the finite time interval [0, t] 
into K subintervals, so that t = KSt, and let 5t — > 0, — > oo, holding t constant. Then 
we have the following representation solution to (1): 

^-«^» 

J] exp {^{t m L - £ m L^ exp (-| \L - l m \ 2 ^ exp p(0) 



K 

X 

m=l 
K 

X 

m=l 



K 

f[ exp exp (-| |L - £ m | 2 ) exp (-|(C^ " ^C 1 ")) (3) 



where H = H+%[L,tf], and the £ m are complex numbers at the discrete moments of 
time labeled by m. We use the notation \L — £ m \ 2 = (Lr — Re£ m ) 2 + (Lj — Im£ TO ) 2 . That 
this is the solution is readily verified by explicit computation [23]. The solution has the 
form of a "measurement process" of the L's, continuous in time, with "feedback" via the 
terms [t m L - l m tf) [24]. The case of many different Lindblad operators, Lj, is readily 
obtained by taking products over j of the appropriate operators in (3) at each moment 
of time. The ordering of the operators at each moment of time is irrelevant in the limit 
8t — > (although the operators at different times are time-ordered, according to increasing 
m). For future reference we write Eq.(3) in terms of a density operator propagator as 
p(t) = K^piO)]. 

Similarly, the solution to the Ito equation has the explicit representation 

m)) = hm n exp I- (2(tf) m L - lH - (tf) m (L) m ) + (L - (L)J6t 
m=l 



xexp\-^H6tj |V(0)). (4) 

where (L) m denotes the expectation of L at time t m = mSt. Eq.(4) expresses an individ- 
ual history \ip£(t)) as an explicit functional of an individual complex gaussian noise £(t). 
However, to make clearer the connection with Eq.(3) (and ultimately with the decoherent 
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histories approach discussed below), it is useful to present an alternative representation, 
obtained by exchanging the Wiener process £ (t) for the stochastic variable £(i) , defined by 
(£(t) - (L))dt = d£*(t) (Ref.[18]). Eq.(4) thus becomes 

lim f5_t\ K/2 



K 

x II exp^^-^)) <*p(-j \L-tmA exp (~H Stj |V(0)> (5) 
m=l 

where ifg is given above after Eq.(3). AMs a (complex) normalization factor, ensuring that 
(ij;(t)\ip(t)) = 1, and appears because the state is now to be thought of as a functional of 
the stochastic process £(i) instead of \ip) = Again the generalization to many 

Lindblad generators is straighforward. That (4), (5) are solutions to Eq.(2) may be verified 
by explicit computation [23]. This explicit representation of the solution clearly indicates 
that the solutions will tend to exhibit localization in the operator L. It also illustrates 
that the solutions have the form of a "trajectory" , or "history" , concentrated around the 
£(t) corresponding to the particular Wiener process 

We may associate a probability with these trajectories. To see this, recall that we 
are working from the outset with a statistical ensemble of solutions to the Ito equation, 
|^), obeying the rules of standard ("classical") probability theory, and the probability 
distribution of the solutions is that implied by the means of d£, etc., given above. Indeed, 
the solution (3) to the master equation, is a sum over £ of \ipg){ipg\, weighted by the 
probability for each solution. As before, it is more useful to express this probability 
distribution as a probability p[£(t)] over the states satisfying Eq.(5). Since Eq.(3) 
must be a mean over \ipi){ipt\, it is easily seen from Eqs.(3) and (5), that the probability 
distribution over the \ipdYs must be p[^(t)J = |A/ji!(£)]| (or what amounts to the same, 
the norm of the state (5) but without the normalization factor N). Note, however, that 
in general (ip£'\ipt) 7^ for £ 7^ £', and thus p[£(t)} may be thought of as a probability for 
histories of values of L only when the solutions \tfjg) are reasonably well-localized in L. 

The Lindblad and Ito equations, (1), (2), are invariant under redefinitions of the Lj's 
by unitary transformations, i.e., Lj — > L 1 - = Ylk=i ^jk^k^ where U^j is a unitary matrix. 
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For the case n = 1, this is just multiplication by a constant phase, and the solutions (3), 
(4), do indeed manifest this invariance [3]. 

Turn now to the decoherent histories approach. The decoherent histories approach 
is a generalization of quantum mechanics to genuinely closed systems, such as the entire 
universe [1,2]. Its aim is to give a predictive formulation of quantum theory applicable 
to closed systems which does not rely on notions of measurement or on the existence 
of an external classical domain. From such a framework, one hopes to understand the 
emergence of the classical world from an underlying quantum one, and the origin of the 
quantum-classical division upon which the Copenhagen interpretation depends. 

In the decoherent histories approach, the mathematical objects one focuses on are the 
probabilities for histories of a closed system. A quantum-mechanical history is defined by 
an initial state p at time t = together with a string of projection operators P ai ■ ■ ■ P Un 
acting at times t\ . . .t n , characterizing the possible alternatives of the system at those 
times. The projections are exhaustive, P a = 1, and exclusive, P a Pf3 = 5 a pP a . Because 
of interference, most sets of histories for a closed system cannot be assigned probabilities. 
The interference between pairs of histories in a set is measured by the so-called decoherence 
functional, 

D(a,g!) = TV [PaM ■ ■ ■ P ai (h) pP^fa) ■ ■ •P a { l (*n)) (6) 

t_ TT-L i_ TTf 

where Pak\Pk) = e h k Pa^ h k , H is the Hamiltonian of the closed system and a denotes 
the string cq • • • a n . When D(a, a') ~ for a ^ a 1 , inteference may be neglected, and one 
may assign the probability p(a) = D(a, a) to the history. Probabilities assigned under this 
condition may be shown to obey the sum rules of probability theory [1] . Sets of histories 
satisfying this condition are called decoherent. Loosely, satisfaction of these conditions 
means that one can "talk about" (i.e., apply classical logic to) the physical properties 
of the system, and think about those properties as if they were definite, without having 
to invoke notions of measurement. Given the Hamiltonian and initial state for a closed 
system, one's initial aim is to determine the strings of projection operators for which the 
decoherence condition is satisfied. 
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The decoherent histories approach is readily applied to the class of open systems con- 
sidered here, i.e., closed systems in which there is a natural separation into a distinguished 
subsystem and the rest. For such systems, a natural set of histories to study are those 
characterized by the properties of the distinguished subsystem at each moment of time, but 
ignoring (i.e., coarse-graining over) the properties of the environment. Histories of this type 
are often decoherent as a result of the interaction between the system and the environment. 
To be precise, consider histories characterized by strings of projections P ai ®I £ ■ ■ ■ P an ®I £ 
at times t\ ■ ■ -t n , where I denotes the identity on the environment. Now assuming that 
the initial density operator factorizes, the trace over the environment may be carried out 
explicitly in the decoherence functional (6), and, in the regime in which a Markovian 
approximation holds, it then has the form 

D(a,a!) = Tr (p^K^JP^ ■ ■ -^[^^b^P^] • • -P^ j) (7) 

Here, K^ +1 is the reduced density operator propagator introduced above, and the trace is 
now over the distinguished subsystem only. 

Given this expression, and given the explicit form of K above, Eq.(3), we may now 
discuss decoherence. For simplicity, consider the case of projections continuous in time in 
the decoherence functional (7). The discrete time version of the decoherence functional 
will contain terms of the form, 

Pa k < fe _J- • '] P a > k = J d 2 i P ak exp (-^\L - i\ 2 ^) [• • •] exp (-||L - £| 2 ) P< (8) 

The operator exp (^—^-\L — £| 2 j is an approximate projection operator in the limit, used 
here, of small 5t. (It is an approximate projector for all 5t if [L,L^] = 0). Now the 
key point is that on the right hand side, we have two different projection operators P ak , 
Pq/ operating on the same Gaussian projection, exp ^— ^\L — £| 2 j . Because Gaussian 
projections are approximately exclusive, the decoherence functional will be approximately 
diagonal in the a^s if we choose the projections P ak also to be Gaussian projections onto 
L: 

P a = exp ^—^K 2 5t\L — K ^5t ^a\ 2 ^j (9) 
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Here, a is a dimensionless and complex continuous label, and fv _1 (5£) _ 2 is the width of 
the projection which will be tuned by the dimensionless parameter k. The approximate 
exclusivity of these approximate projectors means that a, although continuous, has sig- 
nificance only up to order 1. Clearly, with the choice (9), (8) will be very small unless 
i . 

a ~ nSt^i ~ a . Therefore, histories characterized by strings of projections onto L will 
approximately decohere. We thus arrive at our main result: the variables exhibiting local- 
ization in the quantum state diffusion picture are the same as the variables characterizing 
a decoherent set of histories in the decoherent histories approach. 

The operator (9) is an approximate projection operator (under the conditions stated 
above) onto a subset of the spectrum of the (generally non-hermitian) operator L, and the 
label a is complex. Projections of precisely this type have not previously been used in the 
decoherent histories approach, but there is no obvious obstruction to doing so. Indeed, 
the use of such projections is strongly suggested by form of Eq.(8), which arises as a result 
of the invariance of Eq.(3) under L — ■> e*^L. Furthermore, the connection with histories 
characterized by the more familiar types of projectors may be made by specialization to the 
case of a hermitian L (strictly, to the class of operators equivalent to a hermitian operator 
under multiplication by a phase). 

The degrees of localization and decoherence also are related. From the solution to the 
Ito equation (5), the degree of localization is determined by the degree to which L becomes 
concentrated about a particular trajectory £{t). At each time step St, the localization 
width of L is of order St^^. Similarly, from equation (8), one may see that the degree of 
decoherence is also determined by the degree to which L is concentrated about a particular 
value. Loosely, the off-diagonal terms of the decoherence functional are suppressed in 
comparison to the on-diagonal terms (this is the appropriate way to measure approximate 
decoherence [2]) by a factor of order of exp (— k~^\u — c/| 2 x constant), where the constant 
is of order 1. This means that the projectors given by Eq.(9) define an approximately 
decoherent set of histories only if k « 1. 

The degrees of localization and decoherence are related, therefore, in the sense that ap- 
proximate decoherence of histories may be achieved only if the projectors P a characterizing 
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the histories are coarser than the localization width. 

Finally, consider the probabilities for histories. Given approximate decoherence, the 
decoherent histories approach assigns probabilities to histories given by the diagonal ele- 
ments of the decoherence functional (7). Consider Eq.(8), but now with ay. = a',. Under 
the conditions yielding approximate decoherence, the integrand in (8) will be very small 
unless i K~ 1 (St)~'2a. The projection operators P ak then have essentially no effect 
(except to produce a negligible modification of the width of the neighbouring Gaussian 
projectors) so we can drop them. One thus finds that 

P* k K t t l_ i [--}P ak » exp^-^lL-K-^Hajfel 2 ^...] exp (-^\L - K^St^a^ 

(10) 

Using this result, one may then see that the probabilities assigned to these histories in 
the decoherent histories approach have the form of the norm of the state (5) without the 
normalization factor, and with = K^dt^^a^ty). They are therefore exactly the same 
as the probabilities assigned in the quantum state diffusion approach. A more detailed 
account of this work will be presented elsewhere [23] 
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